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Abstract. In view of the numerous examples in the literature it is attempted to 
outline a theory of Heisenberg spin systems possessing dimerized ground states ("DGS 
systems") which comprises all known examples. Whereas classical DGS systems can 
be completely characterized, it was only possible to provide necessary or sufficient 
conditions for the quantum case. First, for all DGS systems the interaction between 
the dimers must be balanced in a certain sense. Moreover, one can identify four special 
classes of DGS systems: (i) Uniform pyramids, (ii) systems close to isolated dimer 
systems, (iii) classical DGS systems, and (iv), in the case of s — 1/2, systems of two 
dimers satisfying four inequalities. Geometrically, the set of all DGS systems may be 
visualized as a convex cone in the linear space of all exchange constants. Hence one can 
generate new examples of DGS systems by positive linear combinations of examples 
from the above four classes. 
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1. Introduction 

Spin systems with exact ground states are rare and hence have found considerable in- 
terest. A trivial case is a system of N unconnected antiferromagnetic (AF) dimers 
which has the product $ of the individual dimer ground states as its unique ground 
state. An interaction between the dimers would in general perturb the ground state, 
but, interestingly, for certain interactions $ remains a ground state. In these cases the 
perturbational corrections of all orders will vanish; the interaction between the dimers 
is, so to speak, frozen for low temperatures. 

Examples of these systems which minimize their energy for a product state $ of 
dimer ground states ("DGS systems") have been constructed and studied in dozens of 
papers. Sometimes DGS systems are also referred to as "valence bond" (VB) systems, 
or, more generally, as "resonating valence bond" (RVB) systems if superpositions of VB 
states are involved. Here I can only mention a small selection of this literature. In the 
seminal papers of Majumdar and Ghosh even s = 1/2 rings are considered with 
constant NN and NNN interactions of relative strength 2 : 1 which possess two different 
DGS according to the shift symmetry of the ring. More precisely, in the second of these 
papers the authors proved as "an interesting by-product" that <3> is an eigenstate of 
the Hamiltonian and conjectured it being a ground state due to numerical studies up 
to 10 spins. One year later, Majumdar [3] mentions a proof of the DGS property for 
Majumdar- Ghosh rings given in a private communication by J. Pasupathy. 

The generalization of these results to arbitrary spin quantum numbers s is due to 
Shastry and Sutherland [I]. A different generalization of Majumdar-Ghosh rings has 
been given by Pimpinelli |3J and re-discovered by Kumar who extended the cou- 
pling within the s — 1/2 ring to 2n-nearest neighbors with strengths J\ = 2n, J2 = 
2n — 1, . . . , J 2n = 1- Though it might not be adequate to call such a spin system still 
a "ring". Already the ring with NNN interaction could better be viewed as a "ladder". 
Ladders with DGS property have also been studied in [Zj and jH]. Other one-dimensional 
models with DGS are certain dimer-plaquette chains, see P and [TP] . 

Another family of two-dimensional DGS systems can be traced back to the work of 
Shastry and Sutherland [TTj on square lattices with alternating diagonal bonds for every 
second square and arbitrary s. These authors also suggest to classify DGS states as a 
"spin liquids" due to their short range correlation. The S(hastry)S(utherland) model is 
physically realized in SrCu2(BOs)2, see jl2] and [T8] . 

Generalizations of the SS model are possible by introducing, additional to the 
nearest neighbor (NN) and diagonal (D) bonds, new types of interaction, namely 
next nearest neighbor (NNN), knight's-move-distance-away (KM) and further-neighbor- 
diagonal (FND), see [13], ^E] and ^Hj. The dimerized state $ is then an eigenstate of 
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the Hamiltonian if the five coupling constants involved satisfy 

Ji ■ J2 : Js : J 4 : J 5 = 1 : 1 : \ : \ : ^ • (1) 

If an inhomogeneous NN coupling is chosen, namely Jo along the dimer bonds and J\ 
for the remaining NN, the choice 

Jo ■ Ji = 7 : 1 (2) 
yields a DGS for the generalized SS model, see [T3j . 

Finally, I mention generalizations of the SS model to arbitrary dimensions ^Zj and 
by superpositions of uniform pyramids (s = 1/2) constructed by Kumar which will 
be reconsidered in section l3~Tl For further related examples, see also [18 . 

In view of the abundance of examples of DGS systems in the literature I am not 
primarily interested in new examples but will try to characterize the class of all ex- 
amples. Unfortunately, I have achieved a complete characterization only in the case of 
classical spin systems. This and the partial results for quantum systems are contained 
in sections 2 and 3. After the general definitions (subsection 2.1) a necessary condition 
is formulated (theorem 1 in subsection 2.2). It says that <£> will be an eigenstate of the 
spin Hamiltonian iff a certain balance condition for the four coupling constants between 
any two dimers is satisfied. For classical spin systems this balance condition can be 
strengthened to the condition of uniform coupling between any pair of dimers (theorem 
2). 

In section 3 I will give sufficient conditions for DGS systems and enumerate four 
special classes of examples. As mentioned above, systems of these classes can be su- 
perposed by positive linear combinations in order to form new DGS systems. Here we 
identify a spin system with its matrix J of exchange parameters and thus understand 
"addition" of systems as the addition of the corresponding matrices. Hence classical spin 
systems and quantum spin systems with any s are not distinguished at the level of J- 
matrices, but, of course, the definition of classical and quantum DGS systems is different. 

Subsection 3.1 describes "uniform pyramids" which are systems of N dimers with 
uniform coupling between all 2N spins (the base of the pyramid) plus an extra dimer 
one spin of which (the vertex of the pyramid) is uniformly coupled to the other 2N spins. 
These pyramids are DGS systems for arbitrary s if the uniform coupling constants are 
suitably chosen. Another class of DGS systems is provided by small neighborhoods 
of unconnected dimer systems (subsection 3.2). The radius r of the neighborhood de- 
creases with s. Although r is not the optimal value, there seems to be a trend that the 
class of DGS systems is shrinking if s increases. This phenomenon can be illustrated by 
examples, see section 4, but has not yet been strictly proven. 
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For small N and s the class of DGS systems can, in principle, be explicitly deter- 
mined. The method is sketched in subsection 13.31 and the result for iV = 2, s = 1/2 is 
given in the form of four inequalities for polynomial functions of the involved four cou- 
pling constants. Recall that for a classical DGS system the coupling between any two 
pairs of dimers must be uniform. Hence it is possible to encode the structure of such a 
system by an N x N matrix G instead of the 2N x 2N matrix J. Then it can be shown 
that the system is a classical DGS system iff this matrix G is positive semi-definite, 

1. e. iff all its eigenvalues (or all its principal minors) are non-negative (theorem 3 in 
subsection 3.4). Moreover, if a coupling matrix J belongs to a classical DGS system, 
then it also belongs to a quantum DGS system for all values of s. Hence theorem 3 
defines a forth class of special DGS systems. 

In section 4 I will consider two examples. The first one (subsection 4.1) consists of 
two dimers which are weakly coupled in a balanced but not uniform way. If e(s) is the 
maximal coupling strength such that the quantum system with spin quantum number 
s is still DGS, then it follows that e(s) — > for s — > oo since the classical system is not 
DGS for all e > 0. 

The second example (subsection 4.2) consists of three dimers and, due to symme- 
try assumptions, normalization and the balance condition, two independent coupling 
constants, say J 2 and J 4 . It is still possible to exactly calculate the convex set of DGS 
systems in the Ji — J4— plane and to illustrate the subsets defined in section 3 for this 
example. 

For the sake of readability of the paper all proofs of the previously formulated 
theorems and propositions are deferred to section 5. In section 6 we investigate 
the geometric structure of the set Cq> of DGS systems, represented as the set of the 
corresponding J-matrices. It is easily shown that C$ is a proper, convex, generating cone 
in the linear space of all symmetric matrices satisfying the balance condition. Moreover, 
we will see that the interior points of C$ are exactly those systems for which $ is the 
unique ground state. Systems at the boundary of C$ have degenerate ground states. 
In particular, the faces of Cq> consist of all DGS systems having the same eigenspace of 
ground states. We close with a summary (section 7). 

2. Definitions and necessary conditions for DGS systems 

2.1. Definitions 

We consider systems of 2N spins with one and the same individual spin quantum number 
s which are grouped into N fixed pairs ("dimers"). To indicate this grouping the spins 
will be denoted by indices fi = (i,S) where i — 1, . . . , N is the dimer index and 5 = 0,1 
distinguishes between the two spins belonging to the same dimer. Further we consider 
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Heisenberg Hamiltonians 

F(J) = ^ J^s, • 8v , (3) 

where s M = (s^, s^, s^) denotes the /z-th spin observable and J the 2N x 2A r -matrix 
of real exchange parameters or coupling constants J^ v . All operators act on a (2s + 1)^- 
dimensional Hilbert space TC = (£) 2 Jl 1 'Hn- If the spin quantum number s is fixed, we 
may identify a spin system with its matrix J. 

Note that the exchange parameters J^ v are not uniquely determined by the 
Hamiltonian H(J) via (JHJ). Different choices of the leading to the same i?(J) will 
be referred to as different "gauges". We will adopt the following gauges: First, the 
antisymmetric part of J does not occur in the Hamiltonian (JHJ). Hence we will follow 
common practice and choose J Mi , = J Ufl . Thus J is a real symmetric matrix . Second, 
since s M -s M = s(s+l) 1 we may choose arbitrary diagonal elements without changing 
HQ), as long as their sum vanishes, Tr J = 0. The usual gauge chosen throughout the 
literature is J m = 0, fj, — 1, . . . , N, which will be called the "zero gauge". In this article, 
however, we will choose the "homogeneous gauge", which is defined by the condition 
that the row sums will be independent of fi, see also |T§] : 

j = J„ = (4) 

V 

Note that the eigenvalues of J may non-trivially depend on the gauge. 

For any dimer with index % let [z'0, zl] denote the ground state of the AF dimer 
Sio ■ Sn which is unique up to a phase and can be written in the form 

[•0, .1] = -=L= J2 i- i r m \m, -m) , (5) 



TCI— — S 

,(3) 



using the eigenbasis \m), m = —s, . . . ,s of . The ground state of a system of 
unconnected AF dimers is the product state 

N 

$ s = (g)[*0,zl], (6) 

i=i 

called the dimerized state; it has the total spin quantum number S = 0. A system J is 
said to admit dimerized ground states (DGS), or to have the DGS property, iff $ s is a 
ground state of H(§), i. e. iff 

($ S |#(J)$ S ) < (tf|ff(J0*) (7) 

for all ^ G TC with 1 | = 1. Let C| denote the set of all spin systems J with the DGS 
property. If the quantum number s is understood, we suppress it and write simply $ 
and C$. 
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Analogous definitions hold for the classical case: Here the spin observables 
are unit vectors, H(J) cl is the Hamiltonian function, defined on the 2iV-fold Cartesian 
product of unit spheres 

N 

V = X S 2 M , (8) 

and $ d C V is the set of all spin configurations satisfying 

$io + s n = o for alH = 1, . . . , N . (9) 

Note that $ d as well as $ s are invariant under rotations. J is said to have the classical 
DGS property iff the minimum of if (J) is assumed for all s G $ c/ . In this case we write 
J G C|f = C|f. 



Necessary conditions for DGS systems 

Whereas a complete characterization of C|, seems to be possible only for small N and 
s or for the classical case s = oo, one can prove a number of partial results, either 
necessary or sufficient conditions for J G C|. 

Necessary conditions of a rather trivial kind can be obtained whenever one finds a 
state $EK such that the r. h. s. of (JJJ) can be explicitely calculated. Less trivial is the 
following 

Theorem 1 $ is an eigenstate of H{S) iff 

JiO,jO + Jiljl — JiO,jl + JilJO (10) 

for all i < j = 2, . . . , N. 

This theorem says that for a DGS system the interaction between any pair of dimers 
has to be balanced in a certain sense: The strength of the inter-dimer parallel bonds 
must equal the strength of the diagonal bonds, but the strength of the dimer bonds may 
be arbitrary. For example, a spin square is never DGS because there are no diagonal 
bonds at all. 

To give another example, is easy to derive from theorem^ the condition ([TJ for $ being 
an eigenstate of the generalized SS model by superposing five suitable dimer pairs of 
the square lattice, see figure [TJ 

If $ is an eigenstate of if (J), it is straight forward to calculate the corresponding 
eigenvalue, since ($|s i( 5 ■ Sj e |$) = for i ^ j, see section . 

Corollary 1 If § is an eigenstate of if (J), then 

N 

#(J)$ = -2 s ( s + l)^J iCM1 $. (11) 

i=i 



Dimerized ground states 



7 




Figure 1. Decomposition of the generalized SS model ^3] ^H] d§| into pairs of dimers 
satisfying I|1U|) . The types of interaction are explained in the Introduction 
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Because (fTU|) is a linear equation, the set of all real, symmetric, traceless 2N x 2N- 
matrices J satisfying (@J) and (fTU|) will be a linear space of dimension ^N(3N — 1), 
denoted by J7$. Further note that the ground state condition ((ZJ) is conserved under 
positive linear combinations of JJ's. Hence the set C$ of DGS systems will form a convex 
cone embedded in the linear space J7$. Further geometrical properties of Cq> will be 
discussed in section |H1 

Another necessary condition for the DGS property is the following: 
Proposition 1 J/J G C$ then Jjo,ii > for all i = 1, . . . , N. 

Note, however, that the coupling between the dimers can be negative and nevertheless 
the system may have a DGS. For example, this may happen for systems close to uncon- 
nected dimer systems, see proposition El 

In the classical case we have similar but stronger results: The balance condition 
fjlUp can be strengthened to a uniform coupling condition: 

Theorem 2 If $ G C% then 

JiO,jO = JiOJl = Jil,jO = Jiljl = e ij (12) 

and 

Jio,a > (13) 
for alii < j = 2,...,N. 

Consequently, we will denote by J7$° = 3t the linear space of all real, symmetric, 
traceless 2N x 2iV-matrices J satisfying (jlj) and (fT2J) . 

3. Sufficient conditions for DGS systems 

We consider four classes of DGS systems which will hopefully cover all examples known 
from the literature by means of positive linear combinations of J-matrices. Note that 
also DGS systems with a different number of dimers can be superposed in this sense. 

3.1. Uniform pyramids 

A uniform pyramid consists of iV dimers with uniform coupling between all 2N spins 
plus another uniform coupling with one extra spin, see figure EJ This extra spin with 
indices (N + 1,0) is considered as a part of another iV + 1-th dimer, in order to make 
it possible to define the dimerized state More precisely, we require 

JiSj e = a > for all z, j = 1, . . . , iV and 5, e G {0, 1} (14) 

except for (z, 5) = (j, e) , 
J iS>N+1 = p > for aU % = 1, . . . , N and 6 G {0, 1} (15) 

Then the following holds: 
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Figure 2. The uniform pyramid is a DGS system if a, {3 satisfy (|16|l . 



Proposition 2 Let J be a uniform pyramid such that 



a I 1 ij s = ~ . 



Then J 6 C$ 



DGS systems generated by uniform pyramids for s = 1/2 have been considered by Ku- 
mar [?]. The s = 1/2 Majumdar-Ghosh ring [2j can be viewed as a superposition of N 
uniform triangles, similarly the SS model [TTj . 

A superposition of uniform pyramids always yields DGS systems with J M „ > 0. But 
this condition is not necessary, as can be seen by the next class of DGS systems. 

3.2. Neighborhoods of unconnected dimer systems 
For matrices J G C<$> we define its spin modulus 

§(J) = liminl (17) 

as the absolute value of the lowest eigenvalue j min of J. It has similar properties as 
a matrix norm. For example, j m m = implies that all eigenvalues of J vanish, since 
Tr J = 0, and hence J = 0. But since, in general, §(— J) ^ §(J), the spin modulus will 
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not be a norm. Nevertheless it can be used to define neighborhoods of matrices J 6 C$ 
because it can easily be shown that 

' <§(!)< II J|| (18) 



2N- V 

holds, where ||J|| = max{ ||Jx|| | ||x|| = 1 } denotes the so-called spectral norm, see [?]. 

° o 

Let J denote the matrix of an unconnected AF dimer system, i. e. Jio,n> for all 

o 

i — 1, . . . , N and all other non-diagonal matrix elements vanish. Of course, JG C$. But 

o 

also a small neighborhood of J still consists of DGS systems. More precisely: 

o 

Proposition 3 Let J be an unconnected dimer system and A = min{Jjo,a|z = 1 . . . iV}. 

o 

Further let J =JJ +A, A G J7$ such that 
Then J e C$. 

This proposition implies that the cone C$ generates J7$, i. e. J<$, = C$ — C$. Hence it is 
not possible to find a smaller subspace of J<& which already contains all DGS systems. 
This stands in contrast to the classical case, see theorem |2j Moreover, the s-dependence 
of the bound in (JT§j) supports the conjecture that the cones C% shrink with increasing s. 



3.3. The case s = 1/2 

In principle, the cone C$ could be exactly determined as follows: Calculate the 
characteristic polynomial p(X) = det(if(JJ) — A 1) of H(S), J G $<$,. Its roots are 
the real eigenvalues E v of H{T}\ one of them, say E , will be the eigenvalue of if (J) 
w. r. t. $. Factor p(X) = q(\)(X — E ). Then J 6 C$ is equivalent to the condition 
that Eq < any root of q(X). One can easily find criteria for this inequality which do not 
assume that the roots of q(X) are known. Consider, for example, the simple case of q 
being quadratic, say q(X) = (A — a) (A — b), a < b. Then obviously E < a iff q(E ) > 
and q'{Eo) < 0. More generally, one can prove the following: 

Lemma 1 Let q{X) = J2f =0 aeX e , ax > 0, have only real roots. Then E < any root of 

{-l) K+n q (n) {E ) > for all n = 0, . . . , K - 1 . (20) 

Thus it is possible to explicitely determine C$ by K = (2s + 1) 2N — 1 inequalities with- 
out calculating the E u . Unfortunately, it is practically impossible to calculate p(X) for 
general J, even by using computer algebra software, except for small values of iV and s. 



I have determined C$ by this method for the simplest case of iV = 2 dimers and 
s = 1/2. The result is the following: 
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Proposition 4 Let N = 2 and s = 1/2 and J G J7$. Rewrite the dinner indices as 

(1,0) = 0, (1,1) = 1, (2,0) = 2, (2,1) = 3. 

JTien J G C<j> i/J irae following four inequalities hold: 

Joi + ^23 - </ 2 - Jl3 > , (21) 

2(Joi + J23) + J02 + J13 > , (22) 

2</oi^23 — (Jl2 — J13) J23 ~ 

J0l((J02-Jl2) 2 + (J02 + Jl3)J23-2J^) >0, (23) 

2(^oi + ^12(^02 + J13) + J23) 

— ^02 — 2J^ 2 — Jf 3 — Joi(^02 + J13 + 6J23) — (J02 + ^13)^23 > . (24) 

We note that the inequalities in proposition 0] can be written in a hierarchical order 
which makes it more convenient to produce examples of DGS systems. To this end we 
define 

Definition 1 



Pi = Joi + ^23 


(25) 


P2 = J02 + Jl3 


(26) 


P3 = J 12 ~ ~ 


(27) 


Pi = J02 - ~ 


(28) 


P5 = Joi 


(29) 



Then the inequalities (|2"T|) to (J2lj) are equivalent to the following: 



< pi (30) 
-2pi < p 2 < Pi (31) 



--y / Pi(2pi-p 2 ) <P3,Pa < ^VPi( 2 Pi-P2) (32) 



Y + 2pi - P2 y-P^ - 2 V ( 2p * PlP2 4p s)( 2 P? - P1P2 - A Pl)j <P5< 

Y + 2 p 1 1 -p 2 { 2pm + \ V ( 2 Pi - ^ - 4pl) (2p? - p lP2 - ^ (33) 

Again, it follows by the convex cone property of C$ that a system of iV dimers 
with s — 1/2 has the DGS property if J can be written as a sum of 4 x 4-submatrices 
satisfying the above inequalities. 

For another application of the direct method to a homogeneous ring of 3 dimers see 
example 14.21 



3.4- The classical case 

In the classical case it is possible to completely characterize all DGS systems. Recall 
that €ij = €ji denotes the uniform interaction strength between two dimers according to 
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Figure 3. The system at the 1. h. s. cannot be a DGS system for fixed e > and 
arbitrary s, since its classical limit at the r. h. s. is not DGS. The classical ground state 
is indicated by small arrows. 



(112)1 . For any J G J7| z define aniVx N- matrix G(J) with entries 

G u = for all i = 1, . . . , N , (34) 
Gij = €ij for all i ^ j — 1, . . . , N . (35) 

Then we have the following result: 

Theorem 3 Let J G J§ , then J G C% iff G(J) is positive semi- definite. 

Recall that G(J) > iff the iV principal minors det(Gij)ij=i t ..^ n > for n — 1, . . . , TV. 
Hence for classical spin systems the DGS property can be checked by testing TV inequal- 
ities. 

This result is also relevant for quantum spin systems since we have the following: 
Proposition 5 C Cf, for all s — |, 1, |, 

4. Examples 

4.1. 

One of the simplest potential DGS systems J(e), see figureEl shows an interesting effect: 
For given s and sufficiently small e it is a DGS system by virtue of proposition El But 
if e > is fixed and s increases, it eventually looses the DGS property Otherwise we 
would get a contradiction since J(e) ^ C|f by theorem |21 and the (normalized) ground 
state energy must converge for s — > oo towards its classical value as a consequence of 
the Berezin/Lieb inequality 20j[21j. 

12. 

Another s = 1/2 system for which C$ can be directly calculated by the method sketched 
in section 13.31 consists of 3 dimers with equal coupling strength Ji and three further 
coupling constants J2, J3, J4, see figure |H According to theorem [T] we must have 
J3 — |(J2 + J4) if this system has a DGS. Hence only three independent variables, 
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Figure 4. A system of 3 dimers with special coupling constants J\ = 1, J2, J4 and 
^3 = \{'h + Ji)- 

say, Ji, J 2 and J 4 are left. These define a 3-dimensional subspace of J7$ and a cor- 
responding sub-cone C$ of C$. C$ can be represented by a convex subset /C$ of the 
J2 — J4— plane which is the intersection of the cone C$ and the hyperplane J\ = 1. The 
permutation (12) (34) (56) swaps J2 and J4 leaving invariant. Hence /C$ is symmetric 
w. r. t. reflections at the axis J2 = 



/C$ is bounded by the lines J4 = 1, J2 = 1 and J4 + J2 = —2 and by two curves 
with small but finite curvature which lie symmetric to the axis J2 = J4, see figure El 
Prominent points of /C$ are 

• (1,1) : Uniform pyramid with j3 = 0, 

• (1.0) and (0,1) : Majumdar/ Ghosh rings, 

• (0,0) : Three unconnected dimers. 

The small neighborhood of (0, 0) which belongs to /C$ according to proposition El is also 
shown. We see that the bound in (fTU|) is far from being optimal in this case. Further we 
have displayed the convex subset C /C$ generated by three dimer pairs according to 
proposition H| 

For larger s > \ the figures /CJ would shrink towards their classical limit which turns 
out to be the line segment 

- \ < h = Ja < 1 • (36) 
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Figure 5. The set /C$ of points with coordinates (J2, J4) such that the system shown 
in figure0]is a DGS system. Certain prominent subsets and points are displayed which 
re explained in section POl 

As mentioned before, the boundary of /C$ corresponds to a degeneracy of the ground 
state. It is easy to identify the "rival" ground states for the straight parts of the 
boundary: If J2 = 1 the second Majumdar/Ghosh dimerized state becomes a rival 
ground state, analogously for J 4 = 1. If the system approaches the line J 2 + J4 = —2 
the ferromagnetic ground state | TTTTTT) becomes the rival ground state. The curved 
parts of the boundary of /C$ correspond to continuously varying families of rival ground 
states. 



5. Proofs 

5. 1 . Proof of theorem Q 

We rewrite the Hamiltonian Q in the form 

F(J) = J n» s m • s v (37) 

1 N 

— Jie,jS S ie ' S jS + 2 JiO,il s iO ' s il (38) 

i^j £,5=0 i=l 

= J2Hij, (39) 

i<j 

where the distribution of the terms of the second sum in (}38|) to the terms is 
arbitrary. We have = Cg> Iv 3 ' such that acts on Hij = Tii® Hj and P 3 ') on 
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the remaining factors. Recall that the dimerized state has the form 

JV 

$ = (g)[zO,il], (40) 

i=l 

where [iO, il] denotes the AF dimer ground state ^ in Hi = Hio ® Ha- 

Lemma 2 $ is an eigenstate of H(J) iff [z0,il] <S> [j'0, jl] is an eigenstate of for all 
i<j = 2,...N. 

Proof: The if-part is obvious. To prove the only-if-part we consider an orthonormal 
basis (e^)x=o,...,L in Hi such that = [i0,zl] where L = (2s + l) 2 — 1. The action of 
Hij on [i0,il] ® [j0, jl] may be written as 



L 



^(4° ® = h^ef® ef . (41) 

A,/x=0 

Let {ef <g> ef)* = eg* <g> ef ] <g> . . . ef ® . . . ef <g> . . . <g> e^. Then 



#(J)$ = ^^-(e^e^)* (42) 
= EE^(^f)* = ^ = ® e o' } )* • (43) 

i<j Xfi i<j 

The terms at the 1. h. s. of (J4*3~j) proportional to (e^ £g> ef)* with A, /i = 1, . . . , L cannot 
cancel since they occur only once in the sum Yli<j- But they don't occur at the r. h. s. of 
(f4*3l . hence = for all A, /x = 1, . . . , L. maps the subspace of Hij with total 
spin quantum number Sij = onto itself, hence Ttnf = ft,^ = for all /i, A = 1, . . . , L. 



Thus only h^p may be non-zero, which means that eg ® eg is an eigenstate of iif, 



In view of lemma El we only need to consider the case of N = 2 dimers with indices 
% < j in the remaining part of the proof. We set $ = [iO, zl] <8> [j0,jl] and rewrite the 
indices according to 

(iO) = 0, (il) = 1, (JO) = 2, (jl) = 3 , (44) 

Consider a modified Hamiltonian of the form 

H'= J^-s v . (45) 

(A«/)e{02,03,12,13} 

Obviously, $ is an eigenvector of iff it is an eigenvector of H', since the difference 
between and H' consists of two dimer Hamiltonians. We note that for fixed 
/i G {0, 1, 2, 3} the three operators 8$,i = 1, 2, 3 form an "irreducible tensor operator", 
i. e. they span a 3-dimensional irreducible subspace with quantum number 5=1. Here 
and henceforward "irreducible" will always be understood as "irreducible w. r. t. the 
product representation of SU(2) in Hij (or similar spaces)". In order to apply the 
Wigner-Eckhardt theorem (WE), see for example [2*2] . we consider states of the form 

^ = ^01® 1p23 , (46) 
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such that ^01 £ ^01 (resp. ^23 £ ^23) belong to irreducible subspaces of 7^ i (resp. H23) 
characterized by their dimension 2S i + 1 (resp. 26*23 + !)• Recall that WE yields 
"selection rules" of the following kind: The matrix element of a component of a 
tensor operator with representation D between two states belonging to irreducible 
representations T)\ and D 2 is nonzero only if D is contained in the product representation 
of Di and D 2 . In the case of irreducible SU(2) representations which are characterized 
by quantum numbers, say, S, Si and S 2 , the above condition simply reads: \Si — S 2 \ < 
S < S1+S2. The dimer ground states [01] and [23] of course span S 2 = representations. 
Then WE yields: 

Lemma 3 ^ only if 3 01 = S 23 = 1. 

Proof: It will suffice to consider only one term of H', for example Jq 2 Sq^ ® s 2 \ since 
analogous arguments apply for the other terms. We skip the factor J 02 and write 

(Vrf ® 4 } \®) = (V>oi ® ^23 ® 4 i} |[01] ® [23]) (47) 
= (^oi|4 J) |[01])^23|4 l) |[23]). (48) 

The first factor vanishes by WE if S01 7^ 1, the second one if S23 7^ 1. I 

Especially, ($|//"'|$) = and hence, if $ is an eigenvector of H' the corresponding 
eigenvalue can only be zero. 

It is well-known that the irreducible subspaces of TC i are eigenspaces of the permutation 
7T i with eigenvalues (— i) 5 oi+ 2s ; analogously for n 23 . For example, if s = 1/2, the S01 = 
singlet subspace of TCox is spanned by the antisymmetric state ^75 (t I — IT), whereas 
the Soi = 1 triplet subspace is symmetric. The terms ■ s v occurring in H' can be 
generated from s ■ s 2 by applying suitable permutations 7r 01 and n 23 . Hence, using 
lemma H3 and the above-mentioned symmetry of i/j under permutations, we obtain 

(m'm= E j^wk-sM (49) 

(M^)e{02,03,12,13} 

= (J02 + J13 ~ J03 ~ Ji 2 )(^|so • s 2 |$) • (50) 

Since (ip\s ■ s 2 \<&) 7^ for a suitable ip we conclude that J 02 + J i3 — J 3 — J12 = iff 
= iff $ is an eigenvector of H' . Together with lemma 121 this concludes the proof 
of theorem^ 

5. 2. Proof of proposition Q] 

Using the same notation as in section 15.11 we conclude by WE that 

([01] ® V23 \sf ® 4?| [01] ® V23) = ([01] |s«| [01])(V 2 3 ^23) = (51) 
for all (fjiv) e {02, 03, 12, 13}. For fixed j G {1, . . . , N} let 

N 

$'= (g)[i0,il]®|s,s) , (52) 

i=l 
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where \s, s) is the ferromagnetic ground state in Ttj. By (Jo"Tj) and analogous equations 
for permuted indices the expectation value ($' |if(J)| contains no interaction terms 
between dimers and thus 



(&\H(3)\&) = 2 



N 



Ji0,iis(s + 1) + JjOjlS 2 



1=1 



> ($|#(J)|$) = 2 



A? 



^ JiO,qs(g + 1) 



i=l 



(53) 



(54) 



In (J5lj) we used the assumption of proposition that $ is a ground state of H(§). It 
follows that JjojiS 2 > —Jjojis(s + 1) and hence Jjoji > 0, which concludes the proof. 

5. 3. Proof of theorem 

First we want to show that (fT0|) also holds for classical DGS systems. 



(55) 



Lemma 4 Let J e C%. Then 

Jij = JiOJO + Jiljl ~~ JiOJO ~ JilJO = 

for all i < j = 2, . . . , N. 

Proof: By assumption, any classical state s satisfying 
&io + s n — o for alH = 1, . . . , N 

minimizes the energy if (J, s). For such states we may write 
if(J, s) = ^2 J ^ s n ' s » 



N 



2 Jie,jS s ie ' s jS ~ 2 JiO,il 

l<i<j<N e,8=0 i=l 

2 (JiOJO + filj'l — ^iOjl — JilJO) s i0 ' s j0 

l<i<j<N 
N 

-2^J, 



(56) 
(57) 
(58) 



■tO,il 



(59) 



i=l 



The function s i— > if (J, s) is constant for all s e "P satisfying (|56)1 . The above equations 
show that also the function (sj) i=li . . tv ^ -^( s ) = Si<i<j<jv Jij s i ' s j * s constant. We 
write 

H(s) =81-^2 J^ s i + ^' Si ' S J ' ( 60 ) 

2<j<N 2<i<j<N 

Since H(s) is independent of Sx, the second factor in the first scalar product in (JoTJj) must 
vanish: J2 2 <j<N J^j s j = 0- By choosing 8j _L s 2 for all j > 2 we conclude J i2 = 0. This 
concludes the proof of the proposition since the numbering of the dimers is arbitrary. ■ 
Next we want to show f)12j) : 
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(61) 



Proof: It will suffice to show J, 



iO.jl 



Jii,io> since Jjojo = Jn,ji follows by applying the 



permutation (01) and the remaining identity Jjoji 
Consider the state s(a) defined by 



Jiojo by 





1 


1 






f 


— cos a 






/ 









f 


— sin a 


\ 











j s il — 




sin a 








1 








— cos a 






V 





J 




\ 





J 




\ 





> 




\ 





/ 



(62) 



and 



s k0 



(o\ 



V 1 / 



for all k 7^ i, j. 



(63) 



Here a is an arbitrary angle < a < it to be fixed later. Let E(a) denote the energy of 
this state. We conclude 

-E(a) = Jj0,il s i0 ' s il + JjO,jl s jO ' s jl + JiO,jO s i0 ' s j0 + Jil,jl s il ' s jl 

+ JiO,jl s iO ' s jl + Jil,jO s il ' s j0 ~~ X] Jko,ki (64) 

fe=l 

(k/o') 

= - (Jio.ii + Jjoji) cos a + (Ja,jo ~ Jioji) sin a - ^2 Jk0 > kl ( 65 ) 

k=l 

= (Jio.ii + Jjoji) (1 - cos a) + (J aj0 - Jioji) sin a + -J5(0) . (66) 

It is obvious that the state s(a) defined by (J52*j) and is a classical DGS for a = 0, 
hence £7(0) is the ground state energy. If J^a + Jjoji < and Jnjo — Jioji ^ we 
may choose the sign of a such that E(a) < E(0) which is impossible due to the last 
statement. Thus we may assume J^a + Jjoji > 0. Hence the expression has its 
minimum at a value a = «o defined by 

Jilj'O — JiOjl 



tana : 

Ji0.il + JjOj'l 

After some algebra we obtain for the corresponding energy 



^E(a ) = (Jio.ii + Jjoji) ( 1 - \/ 1 + 

which is less than ^E(0) if not Jajo — Jioji = 0. 
It remains to show (|13j): 

Lemma 6 7/jG C|? t/ien 

Ji0.il > 
/or alii = 1, . . . ,N. 



Jil 



jO — JiOjl 



J 



i0,il 



J?o 



ioji . 



(67) 



+ 2^(0) > (68) 



(69) 
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Proof: We rewrite if (JJ, s) using (JoTj) . /ij = Jj 0iil and Si = s i0 + s a> = ^ + S>j for 
i<j = 2,...,N: 

N 

F(J, s) = ■ s a 

+ 2 E 

l<i<j<N 
N 

= j2»i(%- 2 )+ E ''A- 5 , • ( 71 ) 

i=l 1<«^J<JV 

The energy of the DGS will be E = —2 X^li A*i- Assume that for some i e {1, . . . , iV} 
we have /Zj < and consider a state s such that Sj = o for all j ^ i and = 2. 
For this state we obtain if (J, s) = 4//j + f^o which contradicts the assumption that the 
system has a DGS. ■ 
This completes the proof of theorem |21 

5. 4 ■ Proof of theorem [3] 

In order to prove theorem |3] we stick to the notation of the last section and rewrite the 
energy of a classical DGS system (fTTj) in the form 

JV 

E(S) = (S\M\S)-2Y,^- (72) 

i=l 

Here S denotes a 3iV-dimensional vector with components 

Q /o(l) c(2) q(3) q(l) q(2) Q (3K , ? oN 

o — \^o 1 , o 1 , o 1 , . . . , o N , o N , o N j , yio) 

and M = G <8>1 3 , where I3 is the 3x3 unit matrix and G has the components 

G u =lk, i = l,...,N (74) 

Gii = e^, l<i^j<N. (75) 

We have to show that the system is DGS iff M, or, equivalently, G is positive semi- 
definite. This follows immediately from M > iff E(S) > —2 YliLi IH = ^0 for all S 
such that < 2 for alH = 1, . . . , N. 

5. 5. Proof of proposition 0] 

Due to ()14j) and (fTB^) the Hamiltonian of a uniform pyramid can be written as 

if (J) = (a-(3)(S 2 N -2Ns(s+l))+/3 [(S N + s N+lfi ) 2 - (2N + l)s(s + 1)] , (76) 

where 

N 1 

s " = E E s ^ ■ ^ 

1=1 8=0 

The eigenvalues of S% are of the form S(S + 1), where S = 0,1, . . . ,2Ns. Obviously, 
the choice S = minimizes (JTBj) if /3 is small enough. In this case $ is a ground state 
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since it has S = 0. If /3 increases, it will eventually reach a certain value /3q where S = 1 
states have the same energy as $. The other values 5* > 1 can be excluded. For j3 = /3 
and s > 1/2 the coincidence of the energies implies 

(a - (3 )(-2Ns(s + 1)) + (3 Q (s(s + 1) - (2iV + l)s(s + 1)) 

= (a - /3b)(l • 2 - 2iVs(s + 1)) + - l)s) - (2N + l)s(s + 1) , (78) 

hence 

@o = i~~r (79) 
1 + s 

and the system has a DGS for (3 < /3q. The value /3o = a in the case s = 1/2 
follows analogously. The only difference is that in this case the minimal eigenvalue 
of (Sn + sn+i,o) 2 will be s(s + 1) = 3/4 and not (s — l)s as in the case s > 1/2. 

5. 6. Proof of proposition [3] 

o 

We consider the matrix J of a system of N uncoupled AF dimers, i. e. 

o o 

Aj =J i0 ,ii> for alH = 1, . . . , iV and J^v— else. (80) 

Without loss of generality we may assume 

Ai = A = min{Ai|z = 1...N} . (81) 

The other assumptions are as in proposition Let 5 min denote the lowest eigenvalue of 
A, hence §(A) = |<5 m , n |. Consider any such that = 1 and \I/ _L $. The two 

o 

lowest eigenvalues of if (J) are 

TV 

E = - 2s(s + 1) \ and (82) 

i=l 
N 

E 1 = - 2s{s + 1) ^ Ai + A x (2 - 2s{s + 1)) (83) 

i=2 

= E + 2Ai . (84) 
Eq belongs to the non-degenerate eigenvector $, hence 

Ei < (V\H(J)\V) . (85) 
The expectation value of the Hamiltonian H(A) can be estimated as follows 

<tf|tf(A)|tf>= ^A^^l^-^l*) (86) 

> 8 min ^2{^\sl\^>) =5 min Ns(s + l) (87) 

> - 2Ai , (88) 
using (fT9"j) and 5 min < since Tr(A) = 0. Thus 

(*\H(3)\V) = (tt|iJ(J)|tf) + (tf|ff(A)|tf) (89) 
> Ei - 2Ai = £ , (90) 
using (}85|) . (}88j) and (JBlj) . This proves that $ is a ground state of if (J). 
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5. 1. Proof of proposition [3] 

Let J G Ci and ||*|| = 1. Then 

N 

(*|ff(J)|*>=X;G«((*|5?|*>-2)+ X) ^((^-Sil*) (91) 

i=l l<i^j<N 

= TrGW + F , (92) 

where W is an iV x N- matrix with elements 

Wij = (V\Si ■ Sj\V) and F = ($|#(J)|$) . (93) 

It can be easily shown that W > 0. Since also G > by theorem El we may conclude 
that Tr(GW) > 0, (V\H($)\V) > F and hence J G C% for all s = 1/2, 1, 3/2, . . .. 

6. Geometrical structure of C$ 

The defining inequalities of J G C$ are 

($1/7(1)1$) < (^l-ff(J)IV') for all V € H, U\\ = 1 . (94) 

For fixed ip, ip ^ ^, (J94"j) defines a closed half space of the real linear space J7#. Thus C$ 
is an intersection of closed half spaces and hence a closed convex cone. For the notion 
of a cone and related notions, see, for example, [23] . 

Moreover, C$ is a proper cone, i. e. J, —J G C$ implies J = 0. Indeed, J, —J G C$ means 
that (<3>|if (J")|<l>) is simultaneously the lowest and the highest eigenvalue of i?(JJ), which, 
due to TriJ(J) = 0, implies #(J) = and J = 0. 

The set of differences C$ — C$ is a linear subspace of J7$; by virtue of proposition El C$ 
contains interior points and hence C$— C$ = J7$. In other words, J7$ is a generating cone. 



A /ace F is a convex subset of a convex set C such that Aci + (1 — A)c2 G F, c\, c-i G 
C, < A < 1 implies c± E F and C2 G F. In words: If a point of F lies in the interior 
of a segment contained in C, then the endpoints of that segment will also lie in F. A 
proper face is a face different from and C. Special faces are the singletons F = {/}, 
such that / never lies in the interior of a segment contained in C; such / are called 
extremal points of C. The dimension d of a face F is defined as the dimension of the 
afline subspace generated by F. Hence extremal points can be viewed as 0-dimensional 
faces. The above definitions follow, for example, [24J; other authors reserve the notion 
of a "face" to the intersection of C with a supporting hyperplane. By virtue of theorem 
0]both notions coincide for the faces of C$. The boundary of a closed convex set C in 
a finite-dimensional linear space is the union of its faces F, excluding F = C. The set 
of all faces of C will form a lattice w. r. t. the set-theoretic inclusion of faces, such that 
Fi A F 2 = Fi fl F 2 , but Fi V F 2 will be the smallest face containing F\ and F 2 , which in 
general larger than F\ U F 2 . 
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It is possible to more closely characterize the faces of Cq>. To this end we note that 
the definition of C$ does not make use of the special form of $ as a product of dimer 
ground states. Hence may defined for arbitrary normalized states <fi G TC. We will, 
however, always restrict this definition to J e J$ such that C$ will be a closed proper 
convex cone in also for the general case. 

Theorem 4 All faces F of C$ are of the form 

k 

F = pi C* , (95) 

i=0 

where O — 0i £ W> an< ^ {4>ii 4>j) = $ij f or a ^ i, j = 0, . . . , k. 

Conversely, any intersection of the form 4.9,5j) will be a, possibly empty, face of C$ . 

The special case F = C<$, is included in ([95]) with = 0. If J lies at the boundary of C$, 
it is contained in a proper face of C$ and hence there will be at least one further ground 
state different from $. Hence we have the following: 

Corollary 2 J is an interior point of &$, iff $ is a non- degenerate ground state of H(J). 
Proof of theorem |U 

(i) Let F = HiLo^i ^ e given with the properties stated in the theorem. We 
want to show that F is a face. Clearly, F is a convex subset of C<j>. Assume 
J = AJi + (1 — A)J2 G F, Jj G C$, < A < 1. If£; = 0we are done, hence we 
may assume k > 0. Let E{ = denote the ground state energies and 



E[ = (0x|i?(Jj)|0i). By assumption, 

A£( + (1 - X)E' 2 = (^(J)!^) = ($|#(J)|$> (96) 

= AEi + (1 - X)E 2 , (97) 

hence 

A(^-E 1 ) + (l-A)(^-E 2 ) =0. (98) 



In ()98)1 both terms are non-negative since E\ and E2 are ground state energies. 
Hence E[ = E\ and <\>\ is a ground state of H(Si). The same holds for all 
<pi, % = 2, . . . , k. This means that JJi G F, analogously J2 G F, and thus F is 
a face of C$. During this proof we will call faces of this kind standard faces. 

(ii) We will prove the following: 

Lemma 7 Each boundary point of a standard face of C$ is contained in a smaller 
standard face. 

Proof: Let F be a standard face of the form (|93)l and L a line in such that 

F n L = {AJ! + (1 - A)J 2 |0 < A < 1}, Jx ^ J 2 . (99) 

Obviously, Ji and J2 are boundary points, and every boundary point of F can be 
obtained in this way. We consider the 1-parameter analytic family of Hamiltonians 
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H(X) = H(XSi + (1 — A)J 2 )- There exists a complete set of eigenvectors ipj(X) 
of if (A) such that the corresponding eigenvalues Ej(X) are analytical functions in 
some neighborhood of A = 0, see, for example, theorem 7.10.1. We may arrange 
the indices such that the first k + 1 eigenfunctions ipo(X), . . . , ipk(X) are degenerate 
and span the same eigenspace of if (A) as the fa, i = 0, . . . , k for all A G [0, 1]. 
According to the cited theorem they must also be degenerate for A G (— e, 0) for 
some small e > since an analytical function is constant on its total domain of 
definition if it is locally constant. But for A G (— e, 0) the ground states of H(X) 
must be other eigenvectors, say ipk+i(X), . . . ,ipi(X), since AJi + (1 — A)J 2 ^ F. By 
continuity, ■ ■ ■ ,ipi(0) are still ground states of if(0). We thus conclude 

J2 e Hi=o C & if = ^fe+i(0), ...,fa = ipi{0). This proves lemma □ 
(iii) It remains to show that each face F of C$ is a standard face. First consider 
the case where F does not consist of a single extremal point. Let J G F be an 
interior point of F and (0i)j=o,...,fc an orthonormal basis of the eigenspace of if (J) 
corresponding to its lowest eigenvalue E . We may set fa = $. Let Ji G F 
be arbitrary, but Ji 7^ J and consider if (A) = if(AJ + (1 — A)Ji), A G R, 
as well as the affine functions A 1— > (fa\H(X)\fa) , i = 0,...,k. They coincide 
at A = 1 since (fa\H(l)\fa) = (fa\H($)\fa) = E . For some neighborhood 
of A = 1, A J + (1 — A)Ji G C$, since J is an interior point of F. Hence $ 
is a ground state of of H(X) and ($|ff (A)|$) is the minimum of the functions 
A 1 — ► (fa\H(X)\fa), i = 0, . . . ,k. This is impossible unless these functions coincide 
for all A G R. It follows that Ji G f^^C^ = K. Hence F is contained in the 
standard face K. This holds also in the case of F consisting of a single extremal 
point. 

If F is not equal to K it must be part of its boundary since F is a face. If J is 

an interior point of F, as before, it must lie at the boundary of K and, by lemma 

[3 must be contained in a smaller standard face K' = Di=o^i' ^ > k. This is a 

contradiction since the ground states of if (J) are fa, ... , fa, see above. 

Thus F = K, i. e. every face of C$ is a standard face and the proof of theorem 0] is 

complete. 

■ 

It will often be convenient to represent the cone C$ by its intersection /C$ = C$ fl P with 
a suitable hyperplane P, see e. g. the example in section [Ol It follows that the faces of 
C$ are in 1 : 1 correspondence with the faces of the closed convex set /C$, except for the 
vertex J = of the cone C<j>. The boundary of /C$ is partly flat, consisting of faces with 
dimension d > 1, partly curved, consisting of extremal points J of /C$ not contained in 
larger faces except itself. For these extremal points J the ground state of if (J) is 
two-fold degenerate: {J} = flLo^r ^ ^ s arrea dy maximally degenerate in the sense 
that any smaller face Hi=o ^« wm ^ e em Pty- 

Finally we note that throughout this section we never used the special structure of $ as 
a product of AF dimer ground states. Hence our analysis will also apply to the cones 
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0$ defined by other ground states 0, but probably the dimerized ground state cone Cj> 
will show the richest structure compared with other cones C^. 

7. Summary 

In this article the classes C|, of DGS systems have been investigated and completely 
characterized in the classical case s = oo. In the quantum case we calculated the linear 
space J7$ generated by C|, but C| itself could only be explicitely determined for small 
iV and s. In the general "lower bound" C C CJ was constructed as the convex 

hull of three (resp. four) special subsets of C| for arbitrary s (resp. s = 1/2). The 
facial structure of the convex cone C% turned out to be anti-isomorphic to the lattice of 
eigenspaces of ground states of if (J), JJ e C%. 

These results could be used for different purposes: They could help to identify concrete 
spin systems as systems possessing dimerized ground states or to better understand 
known examples of DGS systems. A considerable improvement of the given lower bound 
C C C|, seems to be difficult. On the other hand, this article contains elements of a 
"geometry of multi-dimensional level crossing" which could also be interesting in the 
broader context of quantum phase transitions. It seems worth while and possible to 
extend the results on the geometrical structure of CJ. For example, it remains an open 
problem to prove or disprove the conjecture C% C CJ if s > s'. 
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